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We present a theoretical estimate of the atmospheric Newtonian noise due to fluctuations of at- 
mospheric mass densities generated by acoustic and turbulent phenomena and we determine the 
relevance of such noise in the laser-interferometric detection of gravitational waves. First, we con- 
sider the gravitational coupling of interferometer test-masses to fluctuations of atmospheric density 
due to the propagation of sound waves in a semispace occupied by an ideal fluid delimited by an in- 
finitely rigid plane. We present an analytical expression of the spectrum of acceleration fluctuations 
of the test-masses of the interferometer in terms of the experimentally obtainable spectrum of pres- 
sure fluctuations. Second, we consider the gravitational coupling of interferometer test-masses to 
fluctuations of atmospheric density due to the propagation of sound waves generated in a turbulent 
Lighthill process. We present an analytical expression - in the Fourier space - of the spectrum of 
O I acceleration fluctuations of the test-masses of the interferometer. Finally, we discuss the relevance 

of these noise sources in the detection of gravitational waves by comparing the estimated spectral 
densities of Newtonian atmospheric noises considered here to the expected sensitivity curve of the 
VIRGO detector. 
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PACS numbers: Noise (05.40.Ca), sound waves (47.35.Rs), fluid dynamics (47.10.-g), gravitational waves 
(04.30.-w) 



I. INTRODUCTION 
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It is known that in the weak field approximation of Einstein's General Relativity (GR), the linearized equations of 
, GR are hyperbolic, thus implying propagation of gravitational waves with the speed of light The discovery by 
■ Hulse and Taylor of the binary pulsar PSR 1913 — 16, the measured orbital period of which decreases following the 
predictions of GR, constitutes an indirect observational confirmation of the existence of gravitational waves Q . 

The direct detection of gravitational waves is an important goal of contemporary physics and modern technology 
allows to build ground and space based experiments sufficiently sensitive to detect, in a direct way, gravitational waves 
radiated by astrophysical objects [3[. We emphasize that the detection of gravity waves is also important to verify 
the consistency of alternative theories of gravity that are currently not ruled out from a pure theoretical standpoint 
0, Q . Detection of gravitational waves is performed by measuring the relative displacements of several nearly free 
masses which carry the mirrors defining a Michelson interferometer. The measured quantity is the time-dependent 
difference in the lengths of the two orthogonal arms of the interferometer. In principle, this form of antenna can be 
sensitive down to quite low frequencies. In practice, various noise sources will limit the useful bandpass. Several noise 
processes can generate spurious signals in the antenna, masking the effect induced by the gravitational wave. Some 
of these processes, such as seismic and thermal noise, induce displacements of the mirrors ("displacement noise"), 
while others, such as the noise induced by frequency fluctuations of the laser, affect the phase of the optical rays even 
if a real movement of the mirrors is not present ("phase noise"). In order to define the sensitivity of the antenna 
it is necessary to compare the real signal hew (t) to each relevant fake signal /iNoise (t) ■ This comparison is usually 
expressed in terms of the so-called "linear spectral density" which is defined as the square root of the power spectrum 
of the signal In the case of dimensionless amplitudes, the linear spectral density can be therefore expressed in 
units of \/Hz. Considering for instance the French- Italian interferometric gravitational wave detector VIRGO p, it 
turns out that in the low frequency range (below few tens of Hz) the VIRGO sensitivity is limited by the thermal 
noise of the pendulum suspension. Between few tens of Hz and few hundreds of Hz, the dominant mechanism is 
the thermal noise of the mirrors internal modes. At higher frequencies the VIRGO sensitivity curve is limited by 
shot-noise (noise generated by Poisson statistical fluctuations of the number of photons in the light beam). Another 
important source of noise is the so-called gravity-gradient noise, a noise due to fluctuating Newtonian gravitational 
forces that induce motions in the test masses of an interferometric gravitational-wave detector. Gravity gradients are 
potentially important at the low end of the interferometer frequency range, / < 20Hz. Another noise source that is 
important at these frequencies is the vibrational seismic noise, in which the ground's ambient motions, filtered through 
the detector's vibration isolation system, produce motions of the test masses. It should be possible and practical to 
isolate the test masses from these seismic vibrations down to frequencies as low as / ~ 3Hz, but it does not seem 
practical to achieve significant isolation from the fluctuating gravity gradients. Thus, gravity gradients constitute an 
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ultimate low-frequency noise source; seismic vibrations do not. The Virgo Sensitivity Curve is obtained by summing 
up in an incoherent way the spectral noise densities of all the considered noises (seismic noise, shot noise, thermal 
noise, etc.). This incoherent sum implies that the sensitivity curve is obtained by adding together only quadratic 
terms, averaging to zero the "interference" terms. In other words, it is assumed that different sources of noise are not 
correlated among each other and therefore, 
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where, 
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Gravity gradients were first identified as a potential noise source in interferometric gravitational- wave detectors by 
Weiss p}- The first quantitative analyses of such gravity-gradient noise were performed by Saulson and Spero (jf. 
Improvements of Saulson's works were carried out by Thorne and others [1J]|. Finally, T. Creighton at Caltech has 
revisited the relevance of gravity gradient noise due to atmospheric fluctuations [TTj |. 

In this article, we present a theoretical estimate of the atmospheric Newtonian noise generated by fluctuations of 
atmospheric mass densities generated by acoustic and turbulent phenomena and we determine the relevance of such 
noise in the laser-interferometric detection of gravitational waves (l2j ]. The structure of this paper is as follows: in 
Section II, we briefly review Creighton's work on gravity gradient noise generated by atmospheric mass fluctuations. 
In Section III, we present the general expression for the spectrum of an arbitrary Newtonian noise. In Section IV, 
we consider the gravitational coupling of interferometer test-masses to fluctuations of atmospheric density due to 
the propagation of sound waves in a semispace occupied by an ideal fluid delimited by an infinitely rigid plane. We 
present an analytical expression of the spectrum of acceleration fluctuations S$ (fx, fi\ lo) of the test-masses of the 
interferometer in terms of the experimentally obtainable spectrum of pressure fluctuations S p (fx, fa; w). In Section 
V, we consider the gravitational coupling of interferometer test-masses to fluctuations of atmospheric density due to 
the propagation of sound waves generated in a turbulent Lighthill process. We present an analytical expression - in 

the Fourier space - of the spectrum of acceleration fluctuations Ss (kx, £2; ^ of the test-masses of the interferometer. 

Finally, conclusions and final remarks are presented in Section VI. 



II. GRAVITY GRADIENT NOISE DUE TO ATMOSPHERIC FLUCTUATIONS: CREIGHTON'S WORK 



Mass density fluctuations in the atmosphere caused, for example, by acoustic pressure waves or by temperature 
perturbations induce a stochastic gravitational field that couples directly to the test masses of gravitational interfer- 
ometers and produce noise, the so-called atmospheric Newtonian noise. The first to discuss the effects of such noise 
in a quantitative way was Saulson ||. He considered the effects of background acoustic pressure waves on the one 
hand and the motion of massive bodies in the proximity of the interferometer on the other. Saulson concluded the 
atmospheric Newtonian noise he considered would be insignificant even when using advanced gravitational interfer- 
ometric detectors. Later, T. Creighton [ll| revisited the issue of atmospheric Newtonian noise. Among the possible 
sources of mass density fluctuations in the atmosphere, Creighton considered those that he thought potentially most 
significant. The sources he considered were: 

• Low pressure acoustic waves. 

• Fast moving massive bodies in the proximity of the interferometer. 

• Transient atmospheric shock waves. 

• Massive bodies colliding with the ground or the structures housing the interferometer test masses. 

• Perturbations of atmospheric temperature in the vicinity of the detector. 

The largest small scale atmospheric density perturbations are not caused by pressure waves but rather by temper- 
ature perturbations. When heat is transported through a convective atmospheric layer, convective turbulence mixes 
warm and cold air leading to temperature perturbations on all scales up to the order of millimeters. On the time 
scales of interest, these perturbations can be viewed as essentially "frozen" effects in the air mass, while the pressure 
variations are dispersed quickly via sound waves. Therefore, the fluctuations in air density are induced primarily from 
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temperature perturbations which are generally larger than pressure perturbations by several orders of magnitude. 
Although "frozen" in the air mass, these temperature perturbations can cause density fluctuations that vary quickly 
in time, dp = —p^jr, when the wind transports them in space. Indeed, this is the primary source of noise in optical 
astronomy. We remark that in this context, it is also important to study the noise generated by acoustic perturbations 
(density fluctuations) that are either completely absorbed or reflected by the ground and/or structures in the vicinity 
of the interferometer. 

Even though temperature perturbations transported via the wind are the dominant source of atmospheric density 
fluctuations, they do not produce significant Newtonian noise at very high frequencies. This is due to the long time 
intervals that packets of warm and cold air stay in the vicinity of the interferometer test masses. A possible exception 
exist when the air flow forms vortices around the structure in which the interferometer is housed. Since such flows 
can produce a noise spectrum that has a maximum around the frequencies of the typical vortex circulating in the 
proximity of the test mass, the presence of temperature fluctuations in the atmosphere transported by the wind near 
the detector can give rise to a non-trivial source of density fluctuations. Unexpected variations of pressure caused 
by atmospheric shock waves can be potential sources of transient signals that can be detected by the experimental 
apparatus, provided such phenomena occurs in the proximity of the detector. These "shocks" are relevant because 
they can produce significant pressure variations on time scales smaller than 0.1 sec. This time scale corresponds to 
the smallest value of the band-pass of the majority of current interferometric detectors of gravitational waves. Such 
shocks are essentially transient phenomena that may produce spurious signals in the detector, rather than raise the 
noise threshold. It would be useful to know the signal to noise ratio (SNR) that various shocks could produce. 

Although atmospheric shock waves are potential sources of spurious signals in gravitational wave detectors, they 
are readily treatable using environmental detectors. If such sensors record a pressure variation larger than some 
millibar on time scales of 50 — 100 milliseconds, then we could expect spurious signals with (adimensional) amplitude 
of the order of 10 -22 in the range of frequencies 10 — 20 Hz. An example of atmospheric shocks is "sonic booms" 
generated by supersonic bodies (for example a supersonic airplane flying over the space around the detector); they 
could obscure the detection of gravitational waves. Even though such events are rare, if not entirely non-existent, they 
accentuate the potential seriousness of shock waves with respect to other weaker acoustic sources or from acoustic 
sources originating at greater distances from the experimental apparatus. 

Another potential source of spurious signals in the interferometer is the Newtonian noise caused by the motion of a 
single massive body in the proximity of the detector, or the collision of such a body with the experimental apparatus. 
This last possibility is particularly worrisome, since the deceleration of the body can produce high frequency signals. 
Next we discuss one mechanisms that generate atmospheric Newtonian noise from background acoustic pressure waves. 

Consider a pressure plane wave with frequency /, propagating at the speed of sound characteristic of the medium, 
c s . To estimate the relevance of this noise requires comparison of the spectral density of the acoustic noise signal 
to the sensitivity curve of the interferometric detector. The spectral density of the particular form of noise being 



where the indices i identify the various test masses of the interferometer, r min is the dead air radius about the ith test 
mass, and Sp^ (/) is the acoustic noise spectrum measured outside the building enclosing the ith test mass. Moreover, 



the ratio of heat capacities at constant pressure and temperature of the air at room temperature. In what follows we 
briefly reiterate the important points that led to ([3]). 

1 . One supposes that the relative fluctuation of air pressure is small: < 1. 

2. Only sound waves in proximity of the interferometer are considered. Since the interferometer is sensitive only 
to movement of test-masses parallel to the arms of the detector, the gravitational acceleration produced by 
pressure waves on the test masses is reduced by a factor cos 9. Here 9 is the angle between the direction of 
propagation of the pressure wave and the arm of the interferometer. 

3. The test masses of the interferometer is inside a structure that in principle can be used to eliminate noise 
within a characteristic distance r m i n from the test masses. In order to take this into account, a function C(x) 
is introduced. The function C(x) depends on the shape of the structure, the ways in which it reflects sound 
waves, and on many other factors. 

4. The interferometer is located on the ground, not in an ideal homogenous empty space. It is assumed that the 
waves are almost completely reflected from the ground. 



considered is given by [Til ] : 
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Exploiting these four points, the gravitational acceleration produced by sound waves in the propagation direction 
z is: 

g z (t) = G z 

*pf*+7iV ( 4 ) 



IPf \ c s j \ 4/ 



On the other hand, the gravitational wave signal h (t) in the interferometer is related to the acceleration of one of the 
test masses according to 

d 2 h_g(t) 

In frequency space this relation becomes 

h(f) = ^~^P- (6) 

~ h ^ = ^^e.c{^ysp( f ). (7) 

Assuming the noise is stationary and that the directions and amplitudes of the wave modes are uncorrelated, we have 



Combining (j4|) with (J6]) we obtain 
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where (• • •) indicates an average over the modes of the plane wave that contribute to the noise. Finally, we assume 
that the test masses in the interferometer are at a distance equal to several coherent lengths such that the noise 
can be considered uncorrelated and therefore, can be summed linearly. Finally using ^cos 2 #) = 1/3 we obtain ([3]). 
Creighton [Tl|] concludes that the acoustic background and the temperature fluctuations produce Newtonian noise 
that is below the sensitivity threshold of even the most advanced interferometric detectors even though temperature 
perturbations transported along streamlines of non-laminar flux could produce noise within one order of magnitude 
from the maximum sensitivity region at 10 Hz. 

A more comprehensive study of such phenomena requires the use of more sophisticated models. Additionally, shock 
waves in the atmosphere could produce potentially significant spurious signals for an advanced interferometer. These 
signals could be monitored from acoustic sensors places outside the interferometer structure. In order to avoid the 
non-negligible noise of massive objects transported by the wind and to preserve the assumption of linear additivity 
of noise (hypothesis of uncorrelated noise), it would be necessary to construct barriers to keep such massive bodies 
at a secure distance from the test masses. One of the main purposes of the Virgo project is that of achieving good 
sensitivity at low frequencies (around 4—10 Hz). In this frequency band the thermal and seismic Newtonian noise 
represent the dominant sources of noise. If the thermal noise could be reduced at low frequencies using cryogenic 
techniques or by using high Q materials, the seismic Newtonian noise would represent the sensitivity limit at these 
frequencies 0. 

The objective of this article is to present an analytical estimate of the atmospheric Newtonian noise generated by 
fluctuations of atmospheric mass densities and to judge the relevance of this noise in the detection of gravitational waves 
using laser-interferometric techniques. With the development of very refined and sensitive experimental techniques it 
is possible to study directly fluctuation phenomena in various areas of physics. In our view, fluctuation phenomena are 
related in a natural way to transport or convection phenomena. Diffusion phenomena and irreversible thermodynamics 
are based on results derived from the theory of fluctuations. We consider non-quantum fluctuations (hco « ksT) 
[l3j |. which is the conventional hypothesis employed in fluid mechanics. 

Additionally, we suppose that the viscosity coefficients and thermal conductivity of the fluid (atmosphere) are 
non-dispersive, that is, they are independent of the oscillations of fluctuation to. The phenomena that cause such 
fluctuations of atmospheric mass density are divided into two groups: 

1. Fluctuations generated from Acoustic Phenomena. 

2. Fluctuations generated from Turbulent Phenomena. 



We study theoretical models corresponding to each of these groups and present estimates of the spectral density of 
the atmospheric Newtonian noise for each case. 
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III. CALCULATION OF NEWTONIAN NOISE 

Before studying these physical phenomena, we first calculate some relevant quantities that will enable us to quan- 
titatively describe the shape of the Newtonian noise. The effect of atmospheric density fluctuations 8p(x, t) on each 
one of the four test-masses that are suspended on the suspension tower of the interferometer, can be evaluate using 
Newton's force law F = ma with 



3(x, t) = -G d 3 x'Sp(x', t) 



(9) 



where x is the effective position of the test-mass. The test-masses are arranged so as to form the two arms of the 
interferometer where these arms are oriented along the x and y directions. Furthermore, since only variations of the 
relative length of the light path are detected, we are only interested in the x and y components of the acceleration. 
The relative difference of the arm length is given by h (t) 



h(t) 



l r 
I 



,(i) 



,X 2 > 



,(4) 



,(3) 



(10) 



where L is the length of the arms of the interferometer. Therefore, we can link the temporal second derivative of h (t) 
to the Newtonian acceleration of the masses test, 



d 2 h(t) 
dt 2 



(2) 



,(4) 



,(3) 



(11) 



Since we are interested in the spectral amplitude of the atmospheric Newtonian noise, let us consider the Fourier 
transform of h (t) , expressed formally as 



h(uj)= / dth (t) 



In the Fourier space, (fTTj) becomes 



uj 2 h (w) = - (a { 2 



7iW s( 2 ) 



z( 4 ) s( 3 ) 



— a; — a 



(12) 



(13) 



Due to the stochastic nature of the density fluctuations, we are concerned with the ensemble average Sh (u>) defined 
as 

S h (u) = (h(uj)h* (oj)Y (14) 
Equation (T14"]) represents the spectrum of Newtonian noise. Using (fT5|) , (TH|) becomes 



Sh (w) 



L 2 u 4 



-(2) 



4 3) 



(15) 



We point out that (|15|) will be simplified in a convenient manner in the cases we consider by assuming conditions of 
homogeneity (invariance under translations in a statistical sense) and isotropy (invariance under rotation in a statistical 
sense) of the correlation functions describing density fluctuations and therefore, the fluctuations of acceleration. 



IV. ACOUSTIC PHENOMENA 



It is the interplay between the compressibility and inertia of the fluid that supports the propagation of sound waves 
(the oscillating motion of small amplitudes in an incompressible fluid is called sound waves) in the medium. We work 
with the linear theory of acoustics since we consider perturbations that are negligible in the equations of motion. 
We consider exclusively the compressibility and inertia of the fluid, but no other property of the fluid. We obtain 
the linearized equations of the acoustic theory in their simplest non-trivial form. We neglect the influence on the 
propagation of sound waves from viscosity, heat conduction and inhomogeneity at the boundaries. We consider the 
propagation of waves in an ideal fluid, with reference to so-called " small motions" . That is to say, we consider a 
linearized approximation of the Euler equation (l3j 



dv 1 - 

— + ( v ■ Vw = — Vp. 
at p 



(16) 



6 



This procedure is legitimate when the perturbations in pressure and density (Sp and Sp, respectively) are sufficiently 
small compared to the equilibrium state (po an d p , respectively). We proceed to linearize (|16p together with the 
continuity equation 

dp 

at 

and the adiabatic condition 



+ V • (pv) = (17) 



(s it is the entropy for unit mass). Letting p = pa + Sp, p = p + 5p where p and po are the density and pressure 
of the fluid at equilibrium respectively, Sp and Sp are the perturbations of density and pressure satisfying Sp « p 
and Sp << po. Defining the velocity potential function <fr such that V(/> = v (in the reasonable hypothesis that the 
velocity field is irrotational, Vxt = 0) ; we obtain the sound wave equation 



c 



where c s = yj (dp/dp) s is the speed of sound and can be expressed as c s = 1 / y/px s where \ s ~ P 1 (9p/dp) s is the 
adiabatic compressibility of the medium. It is straightforward to verify that Sp and Sp satisfy the same wave equation 
once we recognize the relationship between the fluctuation of pressure and density of a compressible medium is given 
by Sp = c 2 s Sp. We emphasize that the wave equation introduced for <j>, Sr or Sp are linear approximations. We will not 
consider non- linear wave propagation phenomena. We recall however, that non-linearities may produce "distortion" 
of the initial signal or even the emergence of discontinuities (shock waves). At this juncture we highlight the two 
fundamental assumptions that enable us to obtain the linearized equations: 

1. The necessary condition for applying the linearized equations of motion to describe the propagation of sound 
waves is that the speed of the fluid particles is small compared to the speed of sound, v « c s . 

2. The propagation of sound waves in an ideal fluid is considered adiabatic. 

It is reasonable to assume that the process of compression and rarefaction occurring locally is adiabatic even though 
as a consequence of this thermal inhomogeneities arise in the medium. Such inhomogeneities would lead to hypothesize 
the existence of heat exchange between adjacent regions at different temperature. It turns out however, that the time 
scale over which appreciable heat exchange occurs (locally) is much longer than the time scale over which the causes 
(local pressure and temperature gradients) that generate the thermal exchange are sustained. In other words, we 
assume 

^oscillation ^"conduction 

(20) 

where T osc in at i on represents the period of the sound wave, while T con duction represents the characteristic time scale 
over which a significant heat exchange in the medium can occur. These two time scales can be numerically estimated, 
yielding 

r „ ~ ( ^ )2 - C " ~W 1HZ Ysec (21) 

Conduction — ^ — ~^J2 U \~~f~ J 

where A is the wavelength considered and \ is the thermometric conductivity of the air. Clearly the conduction time 
^conduction is much greater than any oscillation time T OSC iii a tion of interest to us. 

A. Propagation of Acoustic Waves in a Semispace 

As a first problem concerning the generation of atmospheric Newtonian noise, let us consider the gravitational 
coupling of the interferometer test-masses to fluctuations of atmospheric density due to the propagation of sound 
waves in a semispace occupied by a compressible ideal fluid (atmosphere), where the semispace is delimited by an 
infinitely rigid plane (ground). It is necessary to integrate (|19p with the boundary condition 

d(t> (x, y, z = 0,J) = Q 



d z 
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This condition expresses the fact that the velocity of the ideal fluid on the surface of separation is purely tangential. 
The general solution of this problem is given by 



c/)(x, t) = / d 3 kf(k)u^(x, t) 



where are plane waves defined as 



it£ (x, t) = exp i(k ■ x — uit) 



(23) 



(24) 



The function f(k) is defined such that it satisfies / (k x , k y , k z ) — f (k x , k v , — k z ) (to ensure normalization) as well as 
the boundary conditions of the problem. In order to find expressions for the acoustic Newtonian noise it is necessary 
to calculate the correlation function associated with the velocity potential <f> relative to fluctuations of acceleration of 
the test-masses. From fluid mechanics 1131 we find 



In Fourier space (|25p becomes 



Sp (x, t) 



Sp(x, uj) 



Po d(j> (x, t) 
' c? dt ' 



lUJ Poi 



x, UJ 



(25) 



(26) 



Multiplying (|26[) for its complex conjugate and averaging (ensemble average) we obtain, 



Ss p (xi, x 2 ; id) = (Sp (x\, ui) Sp* (x 2 , uj)) 



2 2 

w Po 



0(xi, Ul)4> (x*2, U>) 



a, (x, t) = -G I d 3 x'Sp (x', t) 



Furthermore, since 



we find, proceeding as in (|26[) . 

(S a ).j (xi, x 2 ; u>) = (on (xi, uj) a*j (x 2) u)) 



\\x-xt 



= G / a Xi / a x 2 3 06,5,, 1, , x 2 ; ui) 

\W - Till - Toll 

ll x l x l|l ll x 2 x 2| 



(27) 



(28) 



(29) 



where (x, t) is the i th component of the fluctuation of acceleration exerted on the test-mass located at x. This 
fluctuation of acceleration is generated by fluctuations of atmospheric mass density occupying volume V. For economy 
of notation we define the quantity 



Replacing ([17]) and (J3UJ) in (|29p . we obtain 
(hi (xi, ui) a* (x 2 , uj)) = j— ^ 



IfII II2/II 



(Pxi I d 3 x' 2 ICij (x\ — x*i, x' 2 — x 2 ) \4> (x'i, uj) (f> (x 2 , uj) 



Therefore the main problem is that of evaluating the expression for (xi, x 2 , uj), 



S<(> (xi, x 2 ; uj) = (6(xx, u))(j> (x 2 , uj)) = d ki / d k 2 uj; i (x x , w) ut (x 2 , uj) (f{ki)f*(k 2 ) 



(30) 



(31) 



(32) 



or in other words, that of evaluating the correlation function ^/ (kij f* {^2^ ^- Assuming this correlation function is 

homogenous and isotropic, that is to say, assuming invariance (in a statistical sense) under translations and rotations, 
we have 



f(h)t(h)) = (2n) 3 S^ (k, - fc 2 ) /|/(&) 



(33) 
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Including the constants of normalization in the definition of f(k) and substituting (|33|) into (|32p we arrive at 



<j) [x u uj) <j) {x 2 , uj) > = / d ku^ (xi, uj) (x 2 , uj) ( f(k) 



(34) 



The problem has now been reduced to the calculation of the quantity^ / ^kj 
and using v = \<j>, we get 

V • v = - J d 3 kk 2 f(k)u n (x, t) . 
The continuity equation (jTTJ) expressed in frequency space reads 

Sp{x, to) = ^-f-v. 



From the general expression for / 



(35) 



UJ 



Assuming fixed uj and using ([35]). equation ([36]) becomes 



(36) 



Sp(x, uj) 



i Po L 



d^tf(k)ut (x) 



(37) 



and therefore 



(Sp (fi, uj) Sp* (f 2 , uj)) 



p 2 oj e 



dog / d^ix^u}^) (f{k)r{p) 



(38) 



Using the invariance under rotation around the z axis and the constraint conditions on the amplitudes f(k), it is 
found that (13"81) becomes 



S Sp 0?i, x 2 ; uj) 



2p 2 uj e 



1 

y d^(e, uj)j 



(zi - z 2 ) 



(zi + z 2 ) 



(39) 



where r = \J \x\ — x 2 ) 2 + (yi — y 2 ) 2 , J- ({,, uj) = (\f (£, w)| 2 \ an d Jo ( z ) is the Bessel function of the first kind 24|. 
Assuming x\ = x 2 = (0, 0, z), equation (|39[) simplifies and becomes 



(z, uj) 



1 

| d£F(£,a;) 



1 + cos 



2w£z 



(40) 



Since cfSip = 5 P , we obtain 



5 P (z, w) 



2p 2 uj e 



1 



1 + cos 



2oj£z 



(41) 



From ((IT)) we get 



dz cos (az) <S P (2, w) 



rf^J 7 (£, uj) / ciz cos az 



1 + cos 



2w£z 



(42) 



that is 



y dz cos (az) 5 P (z, uj) = T (J^- 



(43) 
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We therefore conclude, 

+00 



[ — —a, u) ) = — ^— - / dz cos (az)S p (z, uj) . (44) 
V 2u / ttPqUJ" J 


The spectrum of pressure S p (z, ui) can be obtained experimentally by use of an acoustic detector (microphone) 
directed along different directions (the parameter a accounts for this directional variability of the acoustic detector) . 
Hence, by integrating (|4"4"|) in z, the quantity F w) can be readily calculated. Using T {^a, ui) it is possible 

to estimate the correlation function associated with fluctuations of acceleration of the test-masses. Alternatively, it 
is possible to determine any correlation starting from the correlations at the ground level. A simple way to verify 
this consists of applying a Fourier transform to the variables x, y, t in the wave equations and boundary conditions. 
The problem is then reduced to integration of an ordinary differential equation with initial conditions depending on 
the values of pressure at the ground level. Due to a lack of detailed measurements of the correlations of pressure, 
it is necessary to further simplify the model. For example, we can start again from the wave equation for pressure 
fluctuation and Fourier transform only the time variable, 



V + %r)6p(?, w) = 

(m&in) = (45) 

- [Sp (r, w)] z=0 = <5p oxp . (x, y, uj) 

In the boundary conditions appears measures of the fluctuations at ground level, Sp cxp . The solution of problem (|45[) 
can be written as the sum on modes at fixed frequency uj, 

dp (r, ui) = J -^A (fc, w) e^cos (7**) (46) 



where 7^ = \/ ui 2 /c 2 — k 2 and the integration is extended to all values of k = (k x , k y , 0) such that 7^ is real. Then 
the correlations of pressure fluctuations can be written as 

S p (n, r 2 ; ui) = f |^ J ^ (A (h, J) A* (fc 2 , W ) ) e^e^ cos (7^) cos (7^2) (47) 

The simplifying hypothesis is that the correlations between the amplitudes of modes have the following form 

(A (fci, w) A* (k 2 , «) ) = (2tt) 2 A (uj) 5 (2) (hi - fc 2 ) (48) 

that is to say, different modes are completely uncorrelated and depend only on the frequency uj. From such hypothesis 
we can rewrite the correlations of pressure fluctuations as 



2 1 

S p (fx, f 2 ; ui) = i (^j A(w) J dijijJo (^-Ri^j cos (-p-V 1 - V 2 ^j cos (~^T~^~ 



T] 2 (49) 



where R\ 2 — \J (xi — x 2 ) 2 + (y\ — y 2 ) 2 . At this point we consider a local measurement of the spectrum of pressure 
fluctuations given by, 



1 / \ 2 
1 / UJ 



S p (0, 0, uj) = — I- j AH. (50) 
Using (|49[) and ([50)) we obtain the result 
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Finally, for the correlation between accelerations we find 

(£„)<,. (fx, r 2 ; lo) = (0, 0; w) / dV x / d 'r'.X,; (n. - n . r 2 ) 



(52) 

i 

x J d VV J Q (j^R'n) cos ^v^V") cos (^v/T - 



'/ 2 



where J m (z) is the Bessel function of the first kind given by, 



For to = 0, Jq (z) becomes |24| . 



J (z) = - exp(zzcostf)^ = ]T(-l) fe M-f. (54) 
^° fc=0 ( v 

Equation ([52|l allows to obtain an analytical estimate for the correlation between the Newtonian accelerations for any 
pair of points in terms of the spectrum of pressure fluctuations which can be experimentally determined. Unfortunately, 
we do not have direct measurements of S p for the VIRGO detector. As evident from (|52[) and (I54[) . (<S a )y (rij ^2; 
has a non-trivial dependence on w and therefore to compare the effect of this noise with the sensitivity curve of Virgo, 
it is crucial to focus on the proper frequency-band of the noise considered. Then, we compare the square root of the 
strain amplitude h rss (/) of Virgo with the square root of the strain amplitude of the noise considered. In principle, 
we could use our analytical estimate together with experimental values extracted from the literature [l4| and provide 
numerical evidence leading to a numerical estimate of the relevance of the noise considered. 



V. TURBULENT PHENOMENA 



Before considering our specific problem, it is useful to briefly discuss the main characteristics of turbulence. Consider 
a turbulent flow of a incompressible fluid. A turbulent flow is by definition unstable: a small perturbation will in 
general be amplified due to non-linearities appearing in the equations describing the flow. Furthermore, it is evident 
from a great amount of experimental data the turbulent flow of an incompressible fluid is rotational, that is to say, 
Q — V x v 7^ 0, at least in certain regions of the space. The set of equations that defines this physical system consists 
is the Navier- Stokes [l3[ equation 



dv 



■ v = — Vp + f V v 



where v it is the kinematic viscosity of the fluid and the incompressible condition reads 

V-v = 0. 



(55) 



(56) 



The system should be integrated taking account of the chosen initial and boundary conditions. The Navier-Stokes 
equation encodes all we need to know about turbulence. However, it is essential to have data from experimental 
observations in order to properly understand the phenomena since it is a highly non-trivial task to analytically 
integrate the equation due to its inherent non-linearities. Among the parameters that characterize the turbulent flow, 
only the kinematic viscosity v — ^ (77 is the dynamic viscosity) appears in the Navier-Stokes equations. The unknown 
quantities to be determined are v and ^. Moreover the flow depends - through the boundary conditions - on the 
shape and dimensions of the body being inserted into the fluid in order to break the irrotationality of the velocity field 
thereby giving rise to turbulence. Since generally the shape of the body is assumed given, the geometric properties 
are characterized by a typical linear dimension denoted I. Let it be a typical speed of the principal flow of the fluid. 
Then every flow is specified by three parameters: is, u and I. The only adimensional quantity that can be constructed 
from these three parameters is the so-called Reynolds number 1Z [13 |. 



ul 

K = — 



V )v 



^2 
vV v 



(57) 
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The numerator of (|57p represents the transport term (or inertial term), while the denominator is the viscosity term. 
When the Reynolds number is small (1Z <C 1), it is permissible to neglect the inertial forces and therefore, the Navier- 
Stokes equations can be linearized. On the other hand, when 1Z » 1 (large Reynolds numbers), the inertial forces 
dominate those of viscosity. In such situations, instabilities develop that lead to chaotic motion (turbulence). 

In order to clarify the mechanism that leads to the emergence of turbulence, we introduce the concepts of broken and 
restored symmetries. The concept of symmetry is central in the study of transition phenomena and fully developed 
turbulence. Transformation symmetries are represented by either continuous or discreet invariance groups associated 
with a specific dynamical theory. Let G be a transformation group acting on the spatially periodic and non-divergent 
function v (x, t). G is said to be the symmetry group of the Navier-Stokes equations if, for all the v that are solutions 
of the Navier-Stokes equations, and Wg G G, the function gv is also a solution. The known symmetries of the 
Navier-Stokes equations are [la ]: 

1. Spatial translations (i, x, v) a, v) with a € R 3 . 

2. Temporal translations (t, x, v) — ► (t + t, x, v) with r G R. 

3. Galilean transformations (t, x, v) — > ft, x + Vt, v + V^j with V € R 3 . 

4. Parity (t, x, v) — > (t, — a?, —v). 

5. Rotation (t, x, v) -> (t, i?f, i?u) with i? e 50(3, R). 

6. Scaling (t, x, v) -> (A (1 ~ Q) t, Ax, A a «) with A e R+ and a e R. 

It is useful to observe that 5 is valid in the limit where I — ► oo and that 6 is valid for f — > (that is, for large 
Reynolds numbers). Furthermore, observe that the symmetry v — > — u is inconsistent with Navier-Stokes equation 
except when the non-linear term is negligible. Unlike the Euler equation, the Navier-Stokes equation is not invariant 
under temporal inversion. This last fact is a consequence of the emergence of dissipative phenomena. Finally, 
notice that all symmetries except 6 are nothing but macroscopic consequences of fundamental symmetries of Newtons 
equations describing (in classical approximation) the microscopic molecular motion of the fluid. Before describing 
what happens when the Reynolds number increases in the fluid, we define the concept of spontaneous symmetry 
breaking |l6j . 

It appears that macroscopic systems generally have a smaller degree of symmetry at low rather than high temper- 
atures. The manifested symmetry at high temperatures is generally a property of the microscopic Hamiltonian of the 
system. As such, it cannot cease to exist even when the symmetries associated with the Hamiltonian appear to be 
violated. The question is, where did the symmetry go? For example the microscopic Hamiltonian of a ferromagnet is 
invariant under rotation. Lowering the temperature of the system certainly does not change this fact. What actually 
changes is the manner in which the symmetry manifests itself. It is natural to assume that in the most perfect 
manifestation of symmetry the ground state should be invariant under symmetry transformations. In many cases, 
the system has several equivalent ground states that can be mapped into each other via symmetry transformation. 
However, since the system can actually be in one and only one of these ground states, the symmetry appears broken. 
When the ground state of the system does not share the same symmetry of the Hamiltonian, the symmetry is said to 
be spontaneously broken. 

Returning to the issue of turbulence, it is experimentally observed that when the Reynolds number increases, the 
symmetries permitted by the equations of motion and the boundary conditions are subsequently broken. However, 
for very high Reynolds numbers there appears to be a tendency to restore the symmetries (in a statistical sense) far 
away from the boundaries. That is to say, the symmetries are restored on average but not for a single realization of 
the velocity field of the system. Such turbulence is referred to as fully developed turbulence [15| . Fully developed 
turbulence is turbulence that is free to develop without any constraints. The only possible constraints are the 
boundaries, external forces or viscosity. It is observed that the structures of a flow that develops on scales comparable 
to the dimensions over which the fluid evolves cannot be properly defined as "developed" . For this reason, no real fluid 
- even if has a high Reynolds number - can be "fully developed" on large energy scales. On smaller scales however, 
turbulence will be fully developed provided the viscosity does not play a direct role in the dynamics at such scales. 

The turbulent flow at sufficiently high Reynolds numbers is characterized by an extremely irregular and disordered 
temporal variation of the velocity field at each point. There is experimental evidence that the velocity field of a 
turbulent flow is fractal in nature [13, EH ■ Turbulent fluids seem to have fractal velocity fields in the sense that the 
increments of the velocity field are proportional to the power 1/3 of the increment of space. In a turbulent flow at very 
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large Reynolds numbers, the average quadratic increment of the velocity between two points separated by distance I 
scales approximately with the power 2 /3 of this distance [l8[ 

([6v (l)] 2 )~li. (58) 

We will not proceed further with regard to turbulence and will introduce new concepts and theories as required. From 
a mathematical point of view, the central problem of turbulence theory is that of obtaining statistical solutions of the 
Navier-Stokes equations. The standard techniques of fluid mechanics are not sufficiently powerful to study turbulence. 
In the second half of the last century a formal analogy has been found between the theory of turbulence and quantum 
field theory. In both cases a system of interacting fields is non-linear, in principle with an infinite number of degrees 
of freedom. From here follows the similarity of the mathematical apparatus used in both theories. For example, the 
method of Feynman diagrams fl9j used to represent equations (transition amplitudes, etc.) may also be applied to 
turbulence theory (tree diagrams, etc.). Our approach to the study of turbulence is to treat it as an analytic statistical 
theory (Kraichnan-Orszag) that relies on dimensional considerations and similarity [20| . 



A. Incompressible Turbulence: Lighthill Process 

Very weak turbulence can be described by linearized dynamical equations [2lj | . Let us consider the case in which 
turbulence is relatively weak, but not so much so that we are able to neglect the non-linear terms in the dynamical 
equations. In this context, we consider as the main quadratic effect the generation of sound by turbulence in a 
compressible medium. The production of sound due to the self-interaction of turbulent vortices (this is the main 
interaction that cause second order effects) happens only when the compressibility of the medium is taken into 
consideration. As an example concerning the generation of mass density fluctuations due to turbulent phenomena, 
let us consider the wave equation of linear acoustics with a source term due to turbulence [l3| 

*-?&)''™--%i>£s? v > (59) 

Equation (f59|) describes the so-called Lighthill process. This phenomenon consists of the generation of acoustic noise 
due to the presence of turbulent fluid flow. In other words, turbulence (or more accurately, the fluctuations of 
turbulent velocity) generates sound [l3j . We wish to emphasize two points that are implicit in (|59p : 

1. Equation ([59")) has been obtained by assuming we are dealing with compressible turbulence. A fundamental 
difference between compressible and incompressible turbulence is the following: for compressible turbulence, 
variations in the velocity field imply variations in the local mass density of the fluid. Fluctuations in the mass 
density imply local variations of pressure that lead to the emission of sound waves. For incompressible turbulence 
instead, changes in pressure do not produce changes in density. 

2. In the source term of (f59|) are present only the fluctuating parts of the turbulent velocity field ( v = u + v 1 , 
where u is the average velocity field and v' is the fluctuating velocity field). The quantity pv^Vj is called the 
Reynolds strain. The laminar portion of the flow, if one exists, does not play any role in generating turbulent 
fluctuations that give rise to acoustic noise in the Lighthill process. 



For reasons of analytical complexity, we will not solve (|59p . We will instead study the simplified equation 

2 1 3 2 \ r „ , d 2 



where the source term is due to incompressible turbulence, keeping the hypothesis of compressibility in the medium 
where the sound waves propagate. Equation l|60p has been derived by neglecting the effects related to viscosity and 
thermal conductivity. Furthermore, it has been assumed that the incompressible velocity fluctuations v' are small 
compared to the average speed of sound c s . Equation (|60|) can then be viewed as describing the generation of sound 
by turbulence with a small Mach number M = U /c s (where U is the characteristic scale of the velocity of the system) 
and not simply by decaying turbulence (that is, turbulence approaching transition to laminar flow). 

Equation (|60|) leads to a number of important consequences. For example the right hand side of the equation is a 
combination of second order derivatives of the field v(x). This means that in absence of boundary conditions, the 
generation of sound waves from turbulence is equivalent to the radiation emitted from a set of acoustic quadrupoles 
(and not by the usual acoustic dipole sources). Therefore, it follows that if there are no boundaries, at small Mach 
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numbers the turbulence does not represent an efficient source of sound. Thus, we assume the "acoustic noise" is 
generated within a bounded region ("turbulent region") of the fluid in which velocity fluctuations occur. We assume 
that the medium surrounding this volume is at rest (this is a much more extended region than the turbulent region 
and is called the "radiation region""). In order to solve (fBT))) we apply analytic statistical theory of turbulence [2l| . 
To simplify calculations we assume we are dealing with fully developed turbulence that is homogenous, isotropic and 
stationary. 

The turbulence is said to be homogenous if all quantities constructed with a set of n points x\,...,x n (at instants 
ti, ...,t n ) are invariant under any translation of this set. In particular [22j 

{u ai (a?i, h) u an (x n , t n )) = (u ai (a?i + y, t\) u a „ {x n + y, t n )) (61) 

where ( • ••) is the usual ensemble expectation value. The turbulence is said to be stationary if all the average 
quantities involved in the n instants (t\, t n ) are invariant under any temporal translation. In particular [221 ] 

(u ai (xi, h) u an (x n , t n )} = (u ai (fi, ti + t) u an (x n , tn + r)) (62) 

Finally, the homogenous turbulence is said to be isotropic if all the average quantities concerning the set of n points 
xi,...,x n (at instants t\,...,t n ) are invariant under any arbitrary rotation. 

It could be argued that there is no turbulent flows that are homogenous or isotropic at large-scale. Isotropy and 
homogeneity can even be debatable at small scales. Nevertheless, these hypotheses enable us to easily exploit analytic 
statistical theories, thereby enormously simplifying the equations of motion. Such theories are quite powerful in the 
sense that they allow to deal with strong non-linearity when the deviation from the hypothesis of non-Gaussianity is not 
large. The point of view adopted in this article is that these (analytic-statistical) techniques describe in a satisfactory 
manner the dynamics of three-dimensional turbulent flows at small scales. The price of this simplification is the gap 
between the situation studied theoretically and that which can be realized in practice. The homogeneity hypothesis 
implies that turbulence is uniform in space and the concept of stationarity can be described as homogeneity in time. 
Isotropy implies there are no preferred directions in space. There cannot be any average velocity in an isotropic 
field since that would immediately imply a preferred direction. Turbulent isotropic and homogenous fields can be 
relatively simple but they are unphysical. In actual velocity fields, energy arises from some average gradient of 
pressure, temperature or mass, and therefore these fields must be anisotropic. Moreover, these fields will be subject 
to specific boundary conditions that imply they are necessarily inhomogeneous. 

Despite these facts, we will consider in this work fully developed, homogenous, isotropic and stationary turbulence. 
That said, we recast (I60|) in frequency space 



-k 2 + — ^ Sp (jt, uj = p^hkjVtj (kj 



with 

j3 z*Ak-x 



(63) 



v 

where V is the volume occupied by the turbulent fluid and p is the average density. We obtain 

2 

5p(k, w) = Po ° s - hk 3 % (k) (65) 
( uj 2 — c 2 k 2 J 

and thus the noise spectrum S p (k, u) associated to the pressure fluctuation becomes, in the frequency space, 

2 4 

S p (k, to) = (dp (k, u) Sp* (k, J) ) = - — PqCs -^ hkjhkm (Vij (k)Vt m (k)) . (66) 
Finally, the noise spectrum becomes 

p 2 

Sp(k, to) — — ^^kikjkik m Ciji m {k) , (67) 

uj 2 — c 2 k 2 
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where 



Ciji m (k) = (ViMvLtf)) ■ (68) 
Our task now is to compute the quantity Ciji m (k), 

Ciji m (k) = [d 3 x[ d'x'e-^'-^B^^ (of, x>) (69) 



with 

im & *) = («< (x) v 3 (x) Vl (f) v m (a?)) . (70) 

The quantity B^\ m (x, x') appearing in (|70|) is known in statistical fluid mechanics as the tensorial statistical moment 
[2(ij |. In the most general case, the stochastic moments of a stochastic vectorial field represent tensors of order k that 
have the form 

Additional limiting conditions on such stochastic fields lead to new symmetry properties of such tensors. The condi- 
tions of homogeneity and isotropy are of particular interest since, in practice, they are the only symmetries considered. 
However, even without these special conditions, the expression of the statistical moment B^ _ cannot be arbitrary 

since it must satisfy special tensorial transformations. The quantity B^ p (a?i, ...,af„) is a tensor of forth order rep- 
resenting a two-point statistical moment. As a working hypothesis we assume that the fields Vi (x) are Gaussian, 
stochastic velocity fields. Under such hypothesis we obtain (in analogy to what is obtained by applying Wick's 
theorem in quantum field theory fl9j ) 

(Vi (x) Vj (x) vi (x) v m (x)) = {Vi (x) Vj (x)} (vi (x) v m (x)) + (v t (x) vi {x)) (vj (x) v m (x)) + (72) 

+ (vi (x) v m (x 1 )) (Vj (x) vi (x')) . 



From turbulence theory it follows that [231 ] 

3 



and 



Vi (x) Vj (x)) = —Sij I dk£ (k) (73) 



( Vi (x) Vj (x')) = / e"*^-*) ( % - ^ ) (74) 



(2tt) j V k 2 

where £ (k) is the Kolmogorov energy spectrum, that in the interval fco << k << k v can be written as 

£(k) =lC sh~i. (75) 

The interval of validity of (|75|) is characterized by fco ~ 2tt/L and k v — TZ 3 ^ 4 /L = (e/i^ 3 ) 1 ^ 4 where L is the linear 
dimension of the volume of the turbulent fluid and /Co is the Kolmogorov constant. The quantity e represents the 
total energy dissipated due to viscous forces and is given by 

dk2v£ (fc) . (76) 



o 



The quantity £ (fc) in (|75p represents the Kolmogorov energy spectrum in the so-called inertial range. This is permitted 
since we are considering a turbulence problem and such turbulence is characterized by inertial modes (fco << k << fc„) 
while the motion of dissipative modes [k v « k « k' u , k' v is the viscous scale) is always laminar (the function of 
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dissipative modes is to absorb energy from inertial modes and dissipate it). Therefore, using (|72]) . (|73[) . (|74| and (J75J) 
the quantity Ciji m (k) becomes 

Cijimik) = ~5ij5 lm j J dkE (k) j y rf 3 ^ y d?x'e-^ s '- s ) + (77) 

Vfco / V V 

2 

+~ (SuS jm + StnSfl) j d 3 x j d 3 x'e~^( s '- s ) ( j dkE (k) 

V V \k 

with I = |a? — x|. For economy of notation, we define 

41 (k) = f d3 *J d 3 x'e- rk i*- s ) (j dk£ (kU (78) 

V V \k / 

and 

2 



41 ( fc) = / av*(^-1 / ^ (*) . (79) 

v y \/c / 

In terms of ([78]) and ([79]) the quantity Ciji m (k) in (|77|) can be written as 

— * 4: ("2") — * 4: ~* 

C l jl m {k) = y5i 3 5lrnl kokv (k) + - (<5«<Sjm + S imSjl)I^ ' kv (k), (80) 

such that 

fc,¥<W, Jim (fc) - ^fc 4 (41 (k) + zrg^ (fc)) . (si) 

Placing (|8Tj) in ([67|l . the noise spectrum <S p (fc, u>) associated to the pressure fluctuation becomes 

S p (k, u>) = ^plci (41 (A) + 22^ (fc)) ^— 7 . (82) 

la; 2 — c 2 fc 2 ) 

Performing a Fourier transform of (5 a )^- (x*i, £2', u>) in (|29[) in the a; variable, we obtain 

5 s (A,w) = (oi (*, to) a* (k, co)) = (^Pj ^S p (k,uj). (83) 

The quantity Ss(k, to) in (j83|) represents the correlation function associated with the fluctuation of acceleration in the 
Fourier space in terms of the spectrum of pressure described in the frequency space. Substituting (]82[) into ()83|1 we 
find 



S a (k, U ) = I (±np G) 2 (41(A) + 241 (k)) ^ (84) 

cj 2 — c?fc 2 



At this point all that remains is to evaluate the integrals 2^ (A) and T 1 ^ (A) . For (A) we find the explicit functional 
form 



41 (*) = (** 1 " *o S ) (kL) - kL cos (kL)} 2 . (85) 

Concerning I*- 2 ^(fc), we notice it can be recast in the following form, 

(2) ,-\ 4tt (4f L 3 ) /Cj^ , , , ....... 









' f sin a 1 


y dlli sin(AJ) j 


/ da s 




1 J as 




\k l I 
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The integral in a can be expressed exactly in terms of incomplete Gamma functions V (a, z) [24 



T(a, z) 



dcf 



/>oo 

/ t^ 1 exp(-t)dt. 

J Z 



(87) 



The remaining expression in (|86p can be numerically integrated for a given choice of the parameters in play. The 
numerical values of the parameters used in our calculations is given in Table 1. The theoretical value of the Kolmogorov 
constant K.q> is borrowed from Reference [25] . 



Parameter 


Symbol 


Value [MKSA-units] 


linear scale of turbulent region 


L 


150 


Newton's constant 


G 


6.67 x 1CT 11 


kinematic viscosity of air at T — 15°C 


va 


1.8 x 1(T 5 


air density 


Po 


1.3 


Reynolds number 


R 


3200 


Kolmogorov constant 


/C 


9.85 



Tabic 1 



Our numerical estimate leads to the square root of the strain amplitude of the atmospheric Newtonian noise Han n (/) 
given by, 



h-ANN (/) 



del' 



lANN 



(/) 



1 iJs s (k, u ). 



Recall that the sensitivity of the Virgo detector is quantified in terms of the square root of the strain amplitude 
h rss (/). For instance, at / ~ 360Hz the best sensitivity at 50% of efficiency is feo% w 1.1 x 10" 20 /\/Hz 

I J /~360-ffz 

[26| . In our work, we have compared Virgo's h rss (/) to Iiann (/) in the frequency range / 6 [4, 10]. We have, 



hANN (/) 
hrss (,/) 



io- 



J /G[4, 10]Hz 



10 



-2(1 



= 10" 3 < 1. 



(89) 



It turns out that the effect of acoustic noise generated in the Lighthill process considered here is at least three orders 
of magnitude below the sensitivity of the VIRGO interferometer. 



VI. CONCLUSIONS 



In this article, we presented a theoretical estimate of the atmospheric Newtonian noise generated by fluctuations of 
atmospheric mass densities due to acoustic and turbulent phenomena and we judge the relevance of such noise in the 
laser-interferometric detection of gravitational waves. First, we considered the gravitational coupling of interferometer 
test-masses to fluctuations of atmospheric density due to the propagation of sound waves in a semispace occupied by 
an ideal fluid delimited by an infinitely rigid plane. We presented an analytical expression of the spectrum of acceler- 
ation fluctuations S a (fx, r* 2 ; CJ ) °f the test-masses of the interferometer in terms of the experimentally determinable 
spectrum of pressure fluctuations S p (r*i, f 2 ; ui). We do not have direct measurements of S p for the VIRGO detector. 
However, values extracted from the literature lead to conclude that the effect would be at least two orders of mag- 
nitude below the sensitivity curve. Second, we considered the gravitational coupling of interferometer test-masses to 
fluctuations of atmospheric density due to the propagation of sound waves generated in a turbulent Lighthill process. 

We presented an analytical expression, in the Fourier space, of the spectrum of acceleration fluctuations S a (ki , fc 2 J to^j 

of the test-masses of the interferometer. We estimated that the acoustic noise generated in the Lighthill process is 
three orders of magnitude below the sensitivity curve of the VIRGO interferometer. 
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